Testing MOND Over a Wide Acceleration Range in X-Ray Ellipticals 
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The gravitational fields of two isolated ellipticals, NGC 720 and NGC 1521, have been recently 
measured to very large galactic radii (~ 100 and ~ 200 kpc), assuming hydrostatic balance of the hot 
gas enshrouding them. They afford, for the first time to my knowledge, testing MOND in ellipticals 
with force and quality that, arguably, approach those of rotation-curve tests in disk galaxies. In 
the context of MOND, it is noteworthy that the measured accelerations span a wide range, from 
more than 10a to about 0.1a o , unprecedented in individual ellipticals. I find that MOND predicts 
correctly the measured dynamical mass runs (apart from a possible minor tension in the inner few 
kpc of NGC 720, which might be due to departure from hydrostatic equilibrium): The predicted 
mass discrepancy increases outward from none, near the center, to ~ 10 at the outermost radii. The 
implications for the MOND-vs-dark-matter controversy go far beyond the simple fact of two more 
galaxies conforming to MOND. 

PACS numbers: 04.50.-h 98. 52. Eh 98.80.-k 



I. INTRODUCTION 

MOND, which posits a departure from standard dy- 
namics as an alternative to dark matter, has been amply 
tested in disk galaxies of all types (for a recent review of 
MOND, see e.g., Ref. [l|). These forceful tests have been 
made possible by the presence of neutral gas, in ordered, 
quasi-circular, rather well-understood motion, observed 
to large radii. Rotation-curve analysis remains the flag- 
ship of MOND testing. 

Testing MOND in individual elliptical galaxies is far 
less advanced. In most ellipticals our only means of prob- 
ing the gravitational field is the study of their stellar dy- 
namics. But, this probe suffers serious uncertainties due 
to our ignorance of the stellar motions (the orbits of stars 
are not circular and can be anisotropic, with unknown 
dependence on radius). For example, Ref. find from 
the stellar kinematics of NGC 7507, some tension with 
the predictions of MOND, but they studied only a rather 
limited set of orbit distributions. In contrast, Sanders 
(2012, shown in Fig. 37 of Ref. [l|) found it easy to get 
a very good match of the same data with the MOND 
predictions, for other orbit distributions. 

Furthermore, this method allows us to probe the field 
only to relatively small radii, i.e., within the stellar com- 
ponent. But ellipticals have, by and large, high mean 
surface densities. Expressed as accelerations, these cor- 
respond to mean accelerations of the order of a or higher, 
within the bulk of the stellar component. MOND then 
predicts that only very modest mass discrepancies should 
be found within these regions. (The term "mass discrep- 
ancy" refers to the ratio of the mass deduced by dynami- 
cal means, assuming Newtonian dynamics— the dynamical 
mass-and the mass detected directly-the baryonic mass.) 
This prediction is indeed born out by the observations. 
But it means that we cannot expect to test MOND, us- 
ing stellar dynamics, in the more interesting regime of 
low accelerations and large predicted mass discrepancies. 

Analysis using strong lensing of quasars avoids the 
first shortcoming but not the second. The "coincidence" 



27ra ps c 2 /d H [3] (d H is the Hubble distance) implies 
that strong lensing cannot occur in the low-acceleration 
regime for a low-redshift elliptical and a source at cosmo- 
logical distances. So, only very mild mass discrepancies 
are predicted by MOND to appear in such studies. Also, 
these studies produce only one, bulk value for the mass 
discrepancy, not a run with radius. This leaves the in- 
terpretation of such analyses at the mercy of knowing 
accurately the stellar mass-to-light ratios of the ellipti- 
cal in question. It is not surprising, then, that one finds 
different groups reaching different conclusions using even 
the same data. For example, Ref. [ij claim meaningful 
(albeit small) tension between the strong-lensing, MOND 
masses, and the stellar masses they deduce in the central, 
high-acceleration regions of several ellipticals, while Ref. 
[5j claim success for the MOND predictions applying the 
same technique to a larger sample, including that of Ref. 
[H (see also Refs. who reached similar conclusions 

as Ref. H). 

In a rare case, the presence of a rotating gas ring 
around the elliptical NGC 2974 permitted [9J to test 
MOND up to a radius of about 12 kpc, using the ring's 
rotation curve. They find good agreement with the pre- 
dictions of MOND for the interpolating function //i (see 
below). But, again, only small mass discrepancies appear 
within the region that is probed. 

Using more luminous test objects, such as planetary 
nebulae, or globular clusters, one can probe individual 
ellipticals to somewhat larger radii. A notable exam- 
ple is the MOND analysis in Ref. [l(| , of data by Ref. 
Ill on planetary nebulae in three ellipticals, reaching 3- 
5 effective radii. But these still do not reach very far, 
where large discrepancies develop, and the method still 
greatly suffers from large uncertainties due to ignorance 
of the orbit distribution. For instance, Ref. [12j used 
the globular-cluster system of NGC 4636 to probe for 
dark-matter and MOND. They find consistency with the 
MOND predictions, but that no strong statement can be 
made since the accelerations probed are not low enough. 

In some ellipticals we can observe x-ray emitting hot 
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gas to relatively large radii. But these tend to be in 
clusters or groups, and, in particular, to be the central 
galaxy of the cluster. The mass discrepancies such anal- 
yses reveal might then be related to those of the cluster 
at large. 

In light of all this, the recent advent of two studies of 
isolated ellipticals, probing their dynamics to large radii, 
is most welcome. The dynamics of NGC 720 [l3| and 
NGC 1521 [lH were studied assuming hydrostatic equi- 
librium in their x-ray emitting gas envelops. The gravi- 
tational field was probed to a little under 100 kpc for the 
former, and a little under 200 kpc for the latter. From the 
deprojected x-ray temperature and density runs, [ill, [HI 
deduced the spherically averaged runs of dynamical mass, 
Md(R), which is equivalent, conceptually, to measuring 
the rotation curve, v(R) = [GM d (R) / 'R] 1 / 2 , or the accel- 
eration field g(R) = v 2 (R)/R. 

The accelerations measured, g, span a range from more 
than 10a o near the center to about 0.1o at the outer 
radii. In this regard, they are as probing as the best of 
rotation-curves. And, while this method brings its own 
systematics (discussed extensively by Refs. [l3l.[l4|-e.g.. 
unaccounted for departures from hydrostatic equilibrium: 
bulk motion of the gas-such as infall, rotation, and tur- 
bulence; uncertainties introduced by having to deproject 
the gas and stellar distributions, etc.) they afford, like 
rotation-curves, analysis free of our ignorance of the test- 
particles orbits. So, all in all, they permit us, for the 
first time, as far as I am aware, to test MOND in indi- 
vidual ellipticals, with almost the same effectiveness as 
in rotation-curve analysis. 1 ' 2 

Here I compare the derived runs of Md with the de- 
tailed MOND predictions for these runs, based on only 
the observed baryon distribution (for reasonable stellar 
M/L values). The analysis and results are described in 
section [ill and discussed in section UTI] 



II. ANALYSIS AND RESULTS 

The dynamical analysis of Refs. [HI, [IH approximates 
the systems as spherical, and so will mine. For a spherical 
system, all existing formulations of MOND as modified 
gravity predict the relation between the measured accel- 
eration, g, and the Newtonian acceleration produced by 
the baryons, g N : 



asymptotic behaviors are /x(x <C 1) ~ x, fi(x 3> 1) ~ 1. 
We can also write 



9n = (J.(g/a )g. 

Here /j,(x) is the MOND interpolating function, aj; 
ing, e.g., in the modified Poisson theory of Ref. [if 



(1) 

sear- 
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1 The deduction of runs of dynamical mass, using gas hydrostatics, 
has, arguably, not yet reached the accuracy of measuring the 
analogous HI rotation curves in spirals. 

2 Weak-lensing (e.g., Ref. [TElO , which does probe large radii, does 
not probe individual galaxies, only some population averages, 
and is still beset by large errors. 



.9 = KW°o)ff* 



(2) 



v(y) is the interpolating function appearing, more di- 
rectly, in QUMOND [17J . The two theories are equivalent 
for spherical systems if the two interpolating functions 
are related by ix(x) — \/v[x(i(x)}. 

The dynamical mass enclosed within radius R is 
M d (R) = G- 1 g(R)R 2 , while g N (R) = GM b (R)/R 2 , 
where M b (R) is the enclosed baryonic mass. Those 
MOND formulations (and the relativistic theories that 
reduce to them in the nonrelativistic limit) thus predict 
that in isolated spherical systems, Md{R) and M b {R) are 
related by 3 



M d (R) 



GMb(R) 



R 2 a 



M b (R). 



(3) 



In Figures [Tj and [H I reproduce from [Tj, [iij . the 
runs of M*(R), the enclosed stellar mass, M g (R), the 
modeled gas mass, and Md{R) deduced for NGC 1521 
and NGC 720, respectively. The stellar contribution is 
for a spherically averaged triaxial model of the depro- 
jected light distribution. The conversion to stellar-mass 
runs assumes their best-fit values of the stellar mass-to- 
light ratios, which they get from a multi-component fit 
to Md(R) (stars + gas + a dark-matter halo): the I- 
band ratio M*/Lj = 2.55(M/L) for NGC 1521 (with 
an adopted distance of 62.1 Mpc), and the -ftf-band ratio 
M*/L K = O.54(M/L) for NGC 720 (adopted distance 
of 25.7 Mpc). For M d (R), [H, give two alternative 
runs: One derived from a model of the x-ray temperature 
and density profiles (extrapolated beyond the measured 
region). The other run of deduced Md they describe as 
based on "a more traditional 'smoothed inversion' ap- 
proach" . 

The MOND predictions of M d (R), from eq.@, are 
shown in the figures, as points, for several values of i?, 
calculated using M b (R) = M*(R) + M g (R). 

I take a = 1.2 • 10~ 8 cm s" 2 (see Ref. [l|). The inter- 
polating function I use is fJ,(x) = fJ-x(x) = x/(l + x), for 
which vx(y) = (1/2) + (1/4 + ljy) 1/2 , which was shown 
in many previous studies 0, [|| Il8rl2"lj | to perform well in 
the range of x values I cover here. The fact that v\ is not 
compatible with solar system constraints for very large y 
(y > 10 6 ) is of no concern here, since I use this form only 
up to y ~ 20 (as in all the above references). For exam- 
ple, using, instead, v{y) = P- L/2 (y) = (1 ~eT y 1 ) _1 (from 
the families defined in Ref. (21|), which is used all along 
in Ref. [1] for MOND rotation curve predictions, would 
have given the same results here: Despite their different 
analytic form, v l/2 and v\ are everywhere equal to better 
than 4.5%, while v 1/2 has the added advantage that it 



In the dark-matter paradigm, the two are not related. 
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FIG. 1: NGC 1521: MOND predictions for the enclosed dy- 
namical masses for several radii (filled squares), shown in com- 
parison with the dynamical masses deduced by Ref . [3] • Ex- 
cept for the MOND points all the rest arc reproduced from 
Ref. [3]: Gray region, with its central solid (black) line is 
the range of dynamical masses based on a model of the tem- 
perature and x-ray density profiles. The open circles are from 
their alternative derivation of the masses based on a "more 
traditional 'smoothed inversion' approach". The dashed (red) 
line is the contribution of the stars M*(R), for their best fit 
AP/Lj = 2.55(M/L) Q , and the dot-dash (blue) line is that of 
the x-ray gas, M g (R). The MOND predictions are given by 
eq.©, with M b (R) = M*(R) + M g (R), and the interpolating 
function V\. 



approaches 1 fast at high y\ so, unlike v\, is consistent 
with solar system constraints. 

The MOND predictions shown in FigJT] use the same 
M* /Li value as adopted by Ref. Q. We see that in 
this case, the MOND predictions agree very well with 
M d {R) deduced by Ref. pj]. For NGC 720, the best 
fit M*/L K = O.54(M/L) found by Ref. Qj] gives a 
MOND-predicted Md{R) that is somewhat lower than 
the deduced run. With the system parameters adopted 
by Ref. [ij], MOND prefers a higher value of M*/L K . 
Without attempting a best-fit procedure (which would 
not be more informative, anyway) I also show in Figj2]the 
MOND predictions for M*/L K = 0.8(M/L) o . For this 
value the agreement is good-as good as that between the 
two alternative runs of Md presented by Ref. [Ril ]. This 
higher value of M* /L K is also in much better agreement 
with values found in earlier, rotation-curve analyses (see 
Fig. 2 of Ref. [Hj], and Fig. 28 of Ref. Q). 



4 More generally, different interpolating functions coincide at high 
and low values of their argument; so we expect some differences 
only in the intermediate-acceleration regime, between about 
10 kpc and a few tens kpc. 



FIG. 2: NGC 720: The same as FigQ](with results from Ref. 
[ll|), for two sets of MOND predictions: one for the same 
M* /L K — O.54(A//L) preferred by the Newtonian analysis 
of Ref. [HI (small rings); the other for M*/L K = 0.8(M/L) o , 
which MOND prefers (filled squares). In the latter case the 
curve for M*(r) should be raised by a factor 1.5. 



III. DISCUSSION 

To better appreciate the import of the results, it is 
helpful to break the full curve comparison into the sev- 
eral salient MOND sub-predictions: (i) In the inner re- 
gions, R < 10 kpc, the measured accelerations are rather 
larger than a [g(10 kpc) « 3a , 2.3a in NGC 1521, 
and NGC 720, respectively, and higher for smaller R]; 
so, MOND predicts only a small mass discrepancy there: 
it predicts that with a reasonable M* /L value the runs 
of Md and Mf, can be made approximately equal in this 
region; which is indeed the case, (ii) MOND predicts 
that, with the preferred M* /L value, the mass discrep- 



ancy should set in where g = GMd(R)/ R 



Indeed, 



this is where the discrepancy becomes ~ 2. (iii) With 
the same M* /L value, MOND predicts correctly the run 
of the mass discrepancy at large radii, where the acceler- 
ations become much smaller than a ; the predicted dis- 
crepancy is Md/Mf, w (oo/gjv) 1 / 2 . Note that a appears 
in (ii) and (iii) in different and independent roles; the 
fact that the same value matches both is nontrivial. 

The two ellipticals considered here differ greatly from 
the extensively studied disk galaxies, in structure, com- 
position (relative contributions of stars, neutral-gas, and 
hot-gas), and formation history. Also, their fields were 
probed by a different method to that of disk galax- 
ies. Hence, the significance of the above results for 
the MOND-dark-matter controversy goes far beyond the 
mere evidence that two more galaxies conform with the 
predictions of MOND. In particular, in the context of the 
dark-matter paradigm it is all the more unexpected that 
the relation between baryons and dark matter, found here 
for the two ellipticals, is described so accurately by the 
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same formula that accounts for disc-galaxy dynamics. 

With all the successes of MOND in accounting for 
the dynamics of galaxies of all types without dark mat- 
ter, we should remember that MOND does not explain 
away completely the mass discrepancy in galaxy clusters, 
where even in MOND some amount of yet undetected 
matter is required (see Ref. [l[ for details). Finding this 
extra mass, or modifying MOND to account for the re- 
maining discrepancy, is a standing challenge for MOND. 

The finer details of the comparison of the MOND pre- 
dictions with the deduced Md runs may be subject to re- 
maining uncertainties (in the distances, the exact value 
of a , systematics, etc.). I want to discuss, specifically, 
the inner regions of NGC 720, where, with the M/L value 
preferred by MOND, the MOND prediction is somewhat 
above the dynamical masses deduced by Ref. |13[. 

In the Newtonian fitting procedure of Ref. [la ], the 
best- fit M* I L value is constrained, as they explain, 
mainly by their deduced Md profile in the inner sev- 
eral kpc. At large radii, Md is dominated by the dark- 
matter halo, which comes with its own adjustable pa- 
rameters. But, it is just these inner regions where the 
x-ray-deduced Md is most susceptible to various system- 
atics (as discussed in Ref. [l3| ). For example, it is there 
that departures from the assumed hydrostatic equilib- 
rium (such as rotation of the gas, turbulence, infall, etc.) 
are most likely. In this connection, the extensive study 
of Ref. [23j], based on a large sample of x-ray ellipticals, 
is most revealing. They find that in the inner parts of el- 
lipticals, where the field is dominated by the stellar mass 
component, there is, generally, little correlation between 
the optical and x-ray morphologies. They conclude that, 
as a rule, "the gas is at least so far out of equilibrium 
that it does not retain any information about the shape 
of the potential, and that X-ray-derived radial mass pro- 
files may be in error by factors of order unity". 5 This 
view is strengthened by the work of Ref. [241 ]. and the 
recent pBT ]. 

In contradistinction, in MOND, M*/L is the only 
"free" parameter that characterizes the mass distribu- 
tion. In the whole region analyzed here, the gas mass 
is subdominant; so, in MOND, M* /L is constrained by 

a Also, departures from spherical symmetry due to the marked 
cllipticity of stellar mass are most important there. 



the whole region of study. So, arguably, more weight 
should be given to the performance of MOND in the re- 
gions beyond several kpc, and the tension in the very 
inner regions of NGC 720 could be easily attributed to 
departures from equilibrium, as found to be ubiquitous 
by Ref. [Hj]. So the MOND M*/L value appears more 
reliable than the value determined by Ref. [l3| . which is 
33 percent lower. In fact, if we permit ourselves even a 
larger discrepancy in the inner few kpc, allowing us an 
even larger M* / L value, based on the inadequacy of the 
x-ray mass determination in the inner region, the agree- 
ment of the MOND predictions with the measure Md 
can be improved in the outer regions. NGC 1521 seems 
to avoid this tension in the inner regions, even assuming 
full equilibrium. 

The accuracy of the MOND predictions of the x-ray- 
deduced fields, over the large range beyond a few kpc in 
NGC 720, and over the whole range in NGC 1521, may 
be taken as evidence that, by and large, the assumption 
of hydrostatic equilibrium is not strongly amiss there. 
Otherwise it would require too much of a conspiracy for 
MOND to correctly predict fictitious Md runs. 

Interestingly, x-ray observations of the inner regions 
of NGC 720 were claimed, in the past, to evince ten- 
sion with MOND: It was claimed in Ref. [2(|, based 
on ROSAT data, that their x-ray isophote ellipticitics 
were too large to be compatible with the MOND field 
of the stellar component alone. Later, it became evi- 
dent that these pronounced ellipticities were largely an 
artifact of point source contamination. Cleaner Chandra 
data [271 showed rather smaller x-ray ellipticities. But 
Ref. [271 ] still maintained that even the remaining ellip- 
ticiies (within ~ 18 kpc), and the degree of misalignment 
they show with respect to the optical isophotes evince the 
presence of a dark matter halo, and-without actually giv- 
ing any calculation in MOND-that a theory like MOND 
cannot account for these. I think that this should be of 
little concern for MOND in light of the finding of Ref. 
[2a] . The case for concluding so is even strengthened by 
what I find in this paper. 

I thank Ofer Yaron for help with the figures. 
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